Absolute magnetometry based on quantum beats in diamond nitrogen-vacancy centers 
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We demonstrate an absolute magnetometer immune to temperature fluctuation and strain inho- 
mogeneity, based on quantum beats in the ground state of nitrogen-vacancy centers in diamond. We 
apply this technique to measure low-frequency magnetic field noise using a single nitrogen- vacancy 
center located within 500 nm of the surface of an isotopically-pure (99.99% 12 C) diamond. The 
photon-shot-noise limited sensitivity achieves 38 nT/vHz for 4.45 s acquisition time, a factor of \/§ 
better than the implementation which uses only two spin levels. For long acquisition times (>10 
s), we realize up to a factor of 15 improvement in magnetic sensitivity, which demonstrates the ro- 
bustness of our technique against thermal drifts. Applying our technique to nitrogen- vacancy center 
ensembles, we eliminate dephasing from longitudinal strain inhomogeneity, resulting in a factor of 
2.3 improvement in sensitivity. 



Negatively-charged nitrogen-vacancy (NV) centers in 
diamond have become an attractive candidate for solid- 
state magnetometry with high sensitivity and nanoscale 
resolution [IH1] , due to their long coherence time [5] and 
near-atomic size. The principle of NV-based magnetom- 
etry is detection of the Zeeman shift of the ground-state 
spin levels. Usually, two spin levels are utilized and the 
presence of a magnetic field induces an a phase shift 
in the spin coherence which can be detected optically 
[Hj. This scheme works well for AC (kHz-MHz) magne- 
tometry [5] and relatively low sensitivity DC field mea- 
surements [H |31 ITHTU] . Sensors based on this technique 
are being developed for applications ranging from neuro- 
science [TUHH], cellular biology [HQS, su P ercon ductivity 
|13| . and nano-scale magnetic resonance imaging |14j . 

Recently, it was discovered that the zero field splitting 
of the NV center ground state is temperature [TS] and 
strain dependent |16j . Consequently, a magnetometer 
using two spin levels is subject to temperature fluctua- 
tion and strain inhomogeneity (if using an NV ensem- 
ble). This limits the magnetometer sensitivity [HI [17] 
(for example temperature fluctuations of 0.01 °C lead to 
fluctuations in the magnetometer reading of ^30 nT) and 
also has implications for quantum information processing 

PS- 

In this Letter, we overcome these issues by exploiting 
the full spin-1 nature of the NV center [2 [TOM] to ob- 
serve quantum beats [351 H3] hi the ground state with a 
beat frequency given only by the external magnetic field 
and fundamental constants. We use a single tone mi- 
crowave field, which transfers all the population into a 
"bright" superposition of the m s = ±1 levels. This tech- 
nique enables measurement of weak magnetic fields at 
the nanometer scale over a broad range of frequencies. 

Quantum beating is a phenomenon of the time evolu- 
tion of a coherent superposition of non-degenerate energy 



eigenstates at a frequency determined by their energy 
splitting. It has wide applications in atomic spectroscopy 
[Ml [35] . and vapor-cell magnetometry [25]. The phe- 
nomenon is closely related to coherent population trap- 
ping, which has been demonstrated in many different sys- 
tems including quantum dots [23 [28] , superconducting 
phase qubits [H] and NV centers [THl I5D]. 

Our quantum-beats magnetometer utilizes a linearly- 
polarized microwave field with frequency / and trans- 
verse amplitude Bmw (perpendicular to NV axis) inter- 
acting with the S = 1 NV ground state (Fig. flji). The 
Hamiltonian describing this interaction is 

JJ/fc = i/ + |l><l|+i/_|-l><-l| (1) 
-fi H ocos(27r/i)(|l)(0| + |0)(1| + |-1)(0| + |0>(-1|), 

where fi^o = SeM-B^Mw/v^ is the undressed Rabi fre- 
quency, [Ib = 13.996 GHz/T is the Bohr magneton, 
g e = 2.003 is the NV electron g— factor, h is Planck's 
constant, and v± is the transition frequency between |0) 
and | ± 1). Here \m s ) denotes the ground state with 
spin projection S z = m s . From Eq. ([IJ, we see that 
the microwave field only drives transitions between |0) 
and a certain superposition of | ± 1), called the bright 
state, \B) = (|1) + |-1))/V2. The orthogonal superposi- 
tion, \D) = (|1) — \-l))/V2, does not interact with the 
microwave field and is therefore called the dark state. If 
Qro ^ |/ - v±\, then Eq. ([TJ describes the Rabi oscil- 
lation between |0) and \B), and the precession between 
\B) and \D) due to the difference of v± can be ignored 
(see Supplementary information). 

Our proposed magnetometer works in the weak field 
and weak transverse strain regime when the transition 
frequencies are [16] 

v ± m D + dh. z ± (g e fi B B z + A\\mj), (2) 
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where D w 2.87 GHz is the ground state zero field split- 
ting, <i" is the axial ground-state electric dipole moment, 
e 2 is the axial electric field (crystal strain), Au = —2.16 
MHz [3T] is the parallel hyperfine coefficient, and mj 
is the spin projection of the 14 N nucleus (I = 1). 
This corresponds to the limit \B\ <C D/g e /j,B ~ 0.1 
T, \g e n B B z + A\\mi\ > {g e ^ B B±_ + A ± mi) 2 /D, and 
\g e HBB z + A\\mi\ 3> |g?j_£jJ, where d± and ej_ are the 
non-axial ground-state electric dipole moment and elec- 
tric field respectively, and A± = —2.7 MHz [32] is the 
perpendicular hyperfine coefficient. Note the last condi- 
tion does not set a minimum detectable magnetic field, 
since for usual diamond samples d±e± is in the kHz range, 
and therefore the condition is always satisfied for at least 
two nuclear sub-levels. 

The experimental scheme, based on Ramsey interfer- 
ometry, is schematically shown in Fig. [JJd. A green laser 
pulse initializes the NV electronic spin into |0). A single- 
tone 7r pulse of sufficient spectral width is then applied 
to transfer the spin into \B). After a free evolution time 
r, the state becomes 

|^(t)) = (e 27IW + T \l) + e 27IW - T \-l))/V2 

= e nl{ "+ +u -^ '{cos[k{u + - V-)t]\B) 

+isin[7r(^ + -^_)r]|D>}. (3) 

We see from Eq. ^ a population evolution between 
\B) and \D) with a beating frequency v+ — V—. Then, a 
second 7r pulse which is phase-coherent with the first tt 
pulse projects the population in \B) back to |0), while the 
population in \D) is trapped. A final green laser pulse 
induces the normalized, ensemble-averaged fluorescence 
signal P(t) oc (1 + J F(r)cos(47r( ffe ^ B B z + A\\mi)r))/2, 

where, for Gaussian decay, F(t) oc e~( T ' with 
the dephasing time. By monitoring P(t) for fixed r w 
(2n+l)/(8(g e fiBB z +A\\mi)) where n is an integer (max- 
imizing the slope of P(r)), we can measure changes in B z . 
Since P(t) only depends on fundamental constants and 
B z , the quantum-beats magnetometer is immune to tem- 
perature fluctuation and strain inhomogeneity, and it is 
absolute without the need of calibration (see Supplemen- 
tary information). 

In comparison, previous magnetometry demonstra- 
tions 2-5 used a large bias magnetic field such that 
coherence between |0) and only one of | ± 1) was se- 
lectively addressed. Broadband magnetometry was re- 
alized using Ramsey interferometry, which begins with a 
green laser pulse used to initialize the spin into |0), fol- 
lowed by a microwave 7r/2 pulse which creates the state 
(|0) + |l))/\/2- After a free evolution of time t, a second 
7r/2 pulse is applied to project the state to |0), which is 
then read out optically. The resulting fluorescence signal 
is P(t) oc (1 + F(t)cos(27T(5t))/2, where 8 = \v+ - f\. 
Since S depends on both D and e z (Eq. ([2])), the mea- 
surement suffers from the temperature dependence of D 
}15j and inhomogenity in t z if using an NV ensemble. 



Our experiments demonstrate that overcoming these 
constraints is critical for high-sensitivity measurement of 
low-frequency magnetic fields. We used an isotopically- 
purified 12 C sample ([ 12 C] =99.99%) [35] to study the 
temperature sensitivity of our quantum-beats magne- 
tometer. The sample has a 500-nm thick isotopically- 
pure layer with [NV]~ 10 11 cm~ 3 grown on top of a 
naturally-abundant substrate with negligible NV density. 
Isotopically purified diamond samples are particularly 
appealing for quantum information and sensing appli- 
cations due to the long spin dephasing times afforded 
by the nearly spinless carbon lattice [5J [TBI I53H5S] . A 
homebuilt confocal microscope was used in the experi- 
ment. Light from a 532 nm laser (~ 1.2 mW) illuminated 
the sample through an oil-immersion objective with 1.3 
numerical aperture, and the fluorescence was collected, 
spectrally filtered, and detected with an avalanche pho- 
todiodc. Pump and probe durations were 2 and 0.3 /^s, 
respectively. A 25 /jm diameter copper wire was attached 
to the surface of the sample to provide square microwave 
field pulses with a Rabi frequency ~ 20 MHz. 

We performed Ramsey interferometry on NV centers 
using both the typical 2-level scheme ({0, 1} basis) and 
quantum-beats detection scheme ({1,-1} basis). A small 
bias field (<2 G) was applied. The spin coherence time 
T£ varies among NV centers in this sample and we chose 
one with relatively long T% . The measured f° r {0; 1} 
basis and {1,-1} basis is 62(2) /j,s and 30(1) /is, respec- 
tively (Fig. ^). 

Using a thermoelectric element, we varied the temper- 
ature of the diamond sample, and performed P(t) mea- 
surements using both the {0, 1} and {1,-1} bases. The 
results are plotted in Fig. [2^,. We see a clear temperature 
dependence in the shape of Ramsey fringes for the {0, 1} 
basis which is not present in the {1,-1} basis. We fit 
the data with a model containing three hyperfine levels 
P(r) — $Z i=1 AjCOs(27rfc'iT -I- fa) + b, where A i7 Vi and 
fa are the amplitude, frequency and phase of the three 
hyperfine oscillations respectively, and b is a constant. 
We used a global fit in which = Vi : r + Ai/(T), and 
Ai, Vi.T , fa and b are fixed for all the temperatures to 
fit for Ais{T) (Fig. [2fj). For the {0,1} basis, assum- 
ing Av{T) = AD(T), we find dD/dT = -78(4) kHz/°C, 
which is consistent with the previous report [15] . Finally, 
we fixed the delay time of the Ramsey interferometer and 
measured the fluorescence level as the temperature was 
varied. As shown in Fig. [2];, the fluorescence level in the 
{0, 1} basis changed significantly and can be well fitted 
with the parameters obtained from fitting the tempera- 
ture dependence of the Ramsey curves. In comparison, 
the change of fluorescence level in the {1,-1} basis is 
about a factor of 7 smaller and does not appear to be 
correlated with the changes in temperature. 

Another advantage of working in the {1,-1} ba- 
sis is the improvement of the magnetometry sensitivity 
by a factor of y/2. Consider the minimum detectable 
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field, b m i n , of a Ramsey-type magnetometer limited by 
quantum-projection fluctuations. It is determined by 
b m in = AN/(Am s g e iiB\dN/dv\), where N = cos(27tvt) 
is the probability distribution difference in the two levels, 
v is the spin-precession frequency in the rotating frame, 
AN is the projection noise, and Am s is the magnetic 
quantum number difference of the two levels. In both 
schemes we only measure the probability distribution in 
two levels, so we can represent the two-level system as 
a spin- 4 system. Then the single-shot projection noise 
is AN = (Aa z ) = (ip\a 2 \ip) - (ip\a z \tp} 2 , where a z is 
the Pauli matrix and ip is the final state. As (ip\cr z \ip) 
is simply the signal TV, we have AN = |sui(27wt)|. In- 
serting this expression into the definition of b m i n , we find 
b mi n{r) = 2 ^ Am J liBC[T)T where C(r) = e^l T ^ is the 
contrast decay due to Gaussian noise. For multiple mea- 
surements, b m i n can be improved by a factor of WT/t, 
where T is the total measurement time. The quantum- 
projection- noise limited sensitivity for our magnetometer 
is thus defined as 



^7min(^") — bminvT 



1 



2TrAm s g e fj, B e 



~(r/T*)\ 



(4) 



The best sensitivity is achieved at r = Although 
T£ in the {1, —1} basis is half of that in {0, 1} basis, Am s 
is twice as big, so rjmin is improved by a factor of \J2 in 
{1,-1} basis (see Supplementary information). 

In our experiment, due to finite photon collection ef- 
ficiency and imperfect spin-state readout, we can only 
measure the photon-shot-noise limited sensitivity (see 
Supplementary information). However, the \[2 improve- 
ment of sensitivity in the {1,-1} basis still persists since 
the contrast of P(t) curves and photon collection effi- 
ciency are the same for the two bases. For T — 5 s, the 
optimal measured sensitivity for the {1,-1} basis and 
{0, 1} basis is 38(3) and 53(4) nT/v 7 !^, respectively. The 
ratio of the two sensitivities is 1.39(15), which is consis- 
tent with the theoretical value of \/2. In comparison, 
the quantum-projection-noise limited sensitivity calcu- 
lated using Eq. Q is 0.79 and 1.12 nT/\/Hz respectively 
(we used the conversion ^/s o y^/v^Hz-)- 

We used both schemes to measure the real noise in the 
laboratory. Using r = 14.5 /xs (29 /zs) for {1, -1} ({0, 1}) 
basis, we repeated the fluorescence measurement in 1- 
second intervals for 50 minutes, now without any active 
temperature control. As seen from Fig. [3^,, for frequen- 
cies near 1 Hz, the noise floor measured by the {1,-1} 
and {0, 1} bases is 50 and 61 nT/yTIz respectively. If we 
ignore dead time from state preparation and readout, the 
sensitivity is 43 and 56 nT/VHz respectively, consistent 
with the photon-shot-noise limit. For lower frequencies, 
the {0, 1} basis suffers more noise which is presumably 
due to laboratory temperature fluctuations. 

To further elucidate this effect, we analyzed the two- 
sample Allan deviation 37J of the fluorescence data, as 



shown in Fig. ^jp. The Allan deviation in the {0, 1} ba- 
sis at long gate time (Tg ato ) levels off and even begins to 
increase at T gat c ~ 100 s, indicating that averaging the 
signal for a longer period of time no longer improves esti- 
mation of a static magnetic field. In contrast, the Allan 
Deviation continues to decrease for the {1,-1} basis up 
to Tgato > 1000 s, indicating this technique is suitable 
for distinguishing nT-scale static fields by using long in- 
tegration times. 

Finally, we studied the effect of strain inhomogeneity 
on a magnetometer employing an ensemble of NV cen- 
ters. We expect that, from Eq. pi), P(t) measurement 
in the {1,-1} basis is insensitive to strain inhomogene- 
ity; however, there will be inhomogenous broadening, and 
consequently reduction of T£ , in the {0, 1} basis due to 
variations in e z for each NV center in the ensemble. We 
used a sample with [ 12 C]=99.9%, which was implanted 
with 10 10 /cm 2 14 N + at an energy of 20 keV and annealed 
at 875 °C for 2 hours, resulting in an NV density of 
~ 5/fim 2 . The laser spot was defocused to illuminate 
a ^2.5 /xm diameter region and the optical power was 
increased to 30 mW to maintain constant intensity. A 
bias field B z sa 0.28 G was applied along the [100] direc- 
tion, such that NVs with different orientation experience 
the same \B Z \. For the {1,-1} basis, we used microwave 
pulses with enough spectral width to cover all three hy- 
pcrfinc levels. For the {0, 1} basis, we detuned the mi- 
crowave frequency and reduced the power to selectively 
address the mi = — 1 level. The peak corresponding to 
the mi = — 1 level of the Fourier transform of P(t) is 
shown in Fig. [4] for both cases. Gaussian fits revealed a 
full-width-at-half-maximum V = 0.09 and 0.12 MHz for 
the {1, —1} and {0, 1} bases, respectively. This indicates 
a factor of ~2.7 increase in spin linewidth due to inho- 
mogenous broadening in the {0, 1} basis, since for a single 
NV r would be half of that of {1, —1} basis. Accordingly, 
we estimate the longitudinal strain inhomcgencity in the 
detected region is ~ 100 kHz. The result indicates a 
factor of 2.3 improvement of sensitivity for the quantum- 
beats magnetometer according to Eq. Q. Our result 
also sheds light on other NV-ensemble applications such 
as quantum memories |38l 139] and frequency references 



In summary, we have demonstrated a magnetometer 
insensitive to temperature fluctuation and strain inho- 
mogeneity based on quantum-beats in NV centers in di- 
amond. The new method uses a similar pulse sequence 
and does not increase the complexity of the magnetome- 
ter. 

Acknowledgement We thank D. Budker, P. R. Hem- 
mer, K.-M. C. Fu, and T. Ishikawa for contributing valu- 
able ideas during the conception of this experiment. K.F. 
acknowledges the support of S. Fan and Stanford Grad- 
uate Fellowship. 



4 



[26 

[27 

[28 
[29 
[30 
[31 

[32 
[33 
[34 
[35 

[36 

[37 
[38 
[39 
[40 



J. M. Taylor, et al., Nat. Phys. 4, 810 (2008). 

J. R. Maze, et al., Nature London 455, 644 (2008). 

C. L. Degen, Appl. Phys. Lett. 92, 243111 (2008). 

G. Balasubramanian, et al., Nature London 455, 648 

(2008). 

G. Balasubramanian, et al., Nature Mater. 8, 383 (2009). 

A. Gruber, et al., Science 276, 2012 (1997). 

L. Rondin, et al., Appl. Phys. Lett. 100, 153118 (2012). 

P. Maletinsky, et al., Nature Nanotechnology 7, 320 

(2012). 

V. M. Acosta, et al., Appl. Phys. Lett. 97, 174104 (2010). 

L. M. Pham, et al., New J. Phys. 13, 045021 (2011). 

L. T. Hall, et al., Scientific Reports 2, 
doi:10.1038/srep00401 

L. P. McGuinness, et al., Nature Nanotechnology 6, 358 

(2011) . 

L.-S Bouchard, et al., New J. Phys. 13 025017 2012. 
M. S. Grinolds, et al., arxiv: 1209.0203 
V. M. Acosta, et al., Phys. Rev. Lett. 104, 070801 (2010). 
F. Dolde, et al., Nature Phys. 7, 459 (2011). 

D. M. Toyli, et al. Phys. Rev. X 2, 031001 (2012). 
P. C. Maurer, et al. Science 336, 1283 (2012). 

E. Togan, et al., Nature 466, 730 (2010). 

P. Huang, et al., Nat. Commun. 2:570 doi: 
10.1038/ncommsl579 (2011). 

F. Reinhard, et al., Phys. Rev. Lett. 108, 200402 (2012). 
F. G. Major, The Quantum Beat: Principles and Appli- 
cations of Atomic Clocks, Springer (2007). 

M. Auzinsh and R. Ferber, Optical Polarization of 
Molecules, Cambridge University, Cambridge, (1995). 

S. Haroche, High-resolution laser spectroscopy, edited by 
K. Shimoda, Springer, Berlin, (1976), pp. 256-313. 

J. N. Dodd and G. W. Series, Progress in atomic 
spectroscopy, edited by W. Hanle and H. Kleinpoppen, 
Plenum, New York, (1978), vol. 1, pp. 639-677. 

S. J. Seltzer, P. J. Meares, and M. V. Romalis, Phys. 
Rev. A 75, 051407(R) (2007). 

B. Michaelis, C. Emary, and C. W. J. Beenakker, Euro- 
phys. Lett. 73, 677 (2006). 

X. Xu, et al., Nature Physics 4, 692 (2008). 

W. D. Kelly, et al., Phys. Rev. Lett. 104, 163601 (2010). 

C. Santori, et al. Phys. Rev. Lett. 97, 247401 (2006). 

B. Smeltzer, J. Mclntyre, and L. Childress, Phys. Rev. 
A 80, 050302 (2010). 

S. Felton, et al. Phys. Rev. B 79, 075203 (2009). 

T. Ishikawa, et al. Nano Lett. 12, 2083-2087 (2012). 

K. Ohno, et al. Appl. Phys. Lett. 101, 082413 (2012). 

N. Zhao, et al. Nature Nanotechnology 
doi:10.1038/nnano. 2012. 152 

K. D. Jahnke, et al., Appl. Phys. Lett. 101, 012405 

(2012) . 

D. W. Allan, Proceedings of the IEEE 54, 221 (1966). 
Y. Kubo et al., Phys. Rev.Lett. 105, 140502 (2010). 
Zhu et al., Nature 478, 221 (2011). 

J. S. Hodges and D. Englund, arXiv: 1109.3241. 




n probe 

or 





20 40 
T ills) 

FIG. 1: a A single-tone microwave pulse interacts with all 
NV center ground-state sublevels. b Ramsey type mag- 
nctomctry using quantum-beats between \B) and \D). c 
Ramsey fringes using the protocol in b. The fitted decay 
envelope yields = 30(1) /is. Inset: Absolute value of 
the fourier transform of the Ramsey fringes. The three 
peaks correspond to the three hyperfine resonances be- 
tween m„ = ±1 levels. 
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FIG. 3: a Measured noise spectrum using the {0, 1} and 
{1, —1} bases. Dashed lines are the noise floor near 1 Hz. 
b Allan deviation of the noise. 
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FIG. 4: Fourier transform of P(t) in both bases for an 
ensemble of NV centers. Each peak corresponds to the 
■mj = — 1 level. Black curves are Gaussian fits. 
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FIG. 2: a Ramsey fringe of the {0, 1} basis (left panel) 
and {1, —1} basis (right panel) for different temperatures, 
b Zero field splitting D dependence on temperature as 
measured in the {0,1} basis (red). The {1,-1} basis 
is immune to changing D (blue), c Fluorescence level 
dependence on temperature for fixed delay time, 18.475 
[is for the {0, 1} basis (upper panel) and 19.720 /js for the 
{1,-1} basis (lower panel). Red curves are measured 
data. Dashed lines indicates the temperature change. 
Solid lines are the calculated fluorescence level for the 
{0, 1} basis. 
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SUPPLEMENTARY INFORMATION 

Detail analysis of possible limitations of our 
approach 

Effect of imperfect n pulse 

In this section, we analyse the influence of the imper- 
fect 7r pulse on the residual temperature-dependence in 
the {1,-1} basis. If the ir pulse does not have enough 
spectral width, then there is some residual population 
in |0) after the ir pulse, which will evolve in the {0, 1} 
basis and thus suffers from the noise from temperature 
fluctuation. Below, we give an estimation of the residual 
population. The Hamiltonian (Eq. (1)) can be simpli- 
fied under the rotating wave approximation and in the 
rotating frame to be time- independent, 



H/h={v + - f)\l){l\ + {v.- f)\-l){~l\ 

■~1}(0| + |0)(-~1| 



(SI) 



~(|1)<0| + |0)(1| 



where |1) = e 27r ^*|l) and |-1) = e 27ri/ *|-l). If we de- 
fine the bright state and dark state in the rotating frame 
as \B) = ^±hl> and \D) = ^p 1 , Eq. Q can be 
rewritten in the basis {|0), \B), \D)} as 
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H/h = -(v+ - v-){\B){D\ + \D)(B\) (S2) 



(\B)(0\ + \0)(B\) 



--K + ^--2/)|0><0|. 
Writing an ansatz for the state 

|*(t)) = Co(t)]0> + c B {t)\B) + c D (t)\D) (S3) 
and inserting it into the Schrodinger's equation, we get 



dc 

— = iitilro 
at 

dc B 
~dT 



inVt m c B + + v~ - 2/)cq, 

i-kCIrqCq - iTi(v + - v-)c D , 



dc D 
~dt~ 



-«7r(^4 



V-)C B - 



(S4) 
(S5) 
(S6) 



Combining Eq. ( S4 ) and Eq. ( S5 ) , we have 

^ 2c • / . os\d C . 2r>2 2o ( \ 

— — in(v + +i/_-2f) — +tt S2 RO c =7r il R o{v + -v_)c D . 

(S7) 

When Qr ^> \v + — Eq. (S7) can be solved us- 
ing perturbation method by expanding the solution cq = 
Cq ' +CQ 1 "* + • • • , where the superindex k represents the or- 



der of 0(( ) fc ). The zeroth order solution gives the 

Rabi oscillation between |0) and \B) with Rabi frequency 



fl R = V^L + i v + 2/) 2 /4. The first order cor- 

rection gives the residual population in |0) after the im- 



(»+-»-)* 



This 



perfect 7r pulse, which is \c \^T~ )l " 
could be taken as a figure-of-merit to estimate the tem- 
perature dependence due to imperfect 7r pulse in {1,-1} 
basis. As an estimation, in our experiment, we have 
v + - v_ ft 4 MHz, fl R ft 20 MHz, / ft (v + + v_)/2. 

The residual population is thus Icq 1 ^ 2 ft 4 x 10 -4 . 



Magnetometer offset due to transverse strain 

The expression for the transition frequency v± (Eq. 
(2)) is modified when |dj_ej_| <C \g e [i B B z + Aym/I is not 
satisfied (the other two conditions for Eq. (2) to hold is 
generally satisfied due to large D). In this case, 



v± = D + dh z ± sJ{g & ii B B z + A ]]mi y 2 + (^e ± ) 2 (S8) 



D + dh z ± (g e fi B B z + Aumi 



2(g e jjL B B z + A\\mi) 



Typical values for single NV centers in our sample are 
d ± e± ft 10 kHz. Thus even for moderate values of \B Z \ ft 
0.5 G, the shift is less than 100 Hz, corresponding to just 
a few nT. 



Temperature dependence due to transverse strain 

The temperature dependence of the transverse strain 
splitting was previously studied for ensembles in a va- 
riety of high-defect-density samples [Slj . In that work 
the term d ± e± was defined as E and the fractional 
temperature dependence was determined to be = 
1.4(3) x 10~ 4 K -1 . Taking the case of E ft 10 kHz and 
\B Z \ ft 0.5 G, th e temperature dependence of the extra 
(d 1 - 



term in Eq. ( S8 1 



is of order 10" 2 Hz/K. 



2(g c fj, B B z +A^m I 

This is more than 6 orders of magnitude smaller than 
dD/dT and undetectable in our experiments. 



Temperature dependence of g e and An 

To our knowledge there have been no published studies 
of the temperature dependence of the NV g-factor and 
hyperfine coupling coefficients. We note that the NV 
electron g factor g e only differs by a tiny amount (~ 
0.03%) from that of a free electron spin due to the small 
spin-orbit interaction in diamond. As a result, we expect 
it is much less sensitive to lattice expansion. 

Measurements of the NV g-factor, g e have been per- 
formed on ensembles at various temperatures from liq- 
uid helium [S2~l to room temperature [S3] . The reported 
values have been consistent, g e = 2.0029, to within the 
experimental uncertainty of ~ 0.0002 [S3] . As a rough 



2 



approximation on an upper bound of the magnitude of 
the effect, we might assume that g e varies linearly by 
.0004 over the temperature interval from 4 to 295 K. In 
this case we find ^fe s» 10 -6 K _1 . The overall temper- 
ature sensitivity of the magnetometer signal would be 
Am s ^ B \B z \^ sa 4Hz/Kfor \B Z \ 0.5 G and Am s = 2. 
This conservative bound is already more than 4 orders 
of magnitude smaller than dD/dT and could be further 
suppressed by operating at lower field. 

Separate measurements of Au for single NV centers 
have been performed at 8 K, Ay (8 K) = -2.169(7) [34], 
and room temperature Ay (295 K) = -2.162(2) [35]. The 
difference in the measured values is within the experimen- 
tal error. A rough estimate can be obtained by taking 
a linear extrapolation, yielding dAnJdT w 24(25) Hz/K. 
This is more than 3 orders of magnitude smaller than 
dD/dT and is not discernable in our experiments. 



T£ of {1,-1} basis and {0,1} basis 

The coherence in the {0, 1} basis can be characterized 
by the average of the phase factor as i ,i = (e^*- 1 ) sub- 
ject to environment noise [SB], where (p(t) = f {v + — f)dt. 
Then the coherence in the {1,-1} basis can be charac- 
terized by the average of the phase factor as £1,-1 = 
(e 211 ^*)). The phase <j>{t) is correlated to the environ- 
mental noise S (for example, through Eq. (2) for mag- 
netic noise), and thus the {1,-1} basis suffers from twice 
as much noise as the {0, 1} basis. When the noise cor- 
relation time is much longer than the dephasing time 
TJ, e.g. Gaussian noise, then 1/T£ is proportional 
to the noise amplitude $71. As a result, in this case 
T*({0, 1}) = 2T*({1,-1}), and we have a V% improve- 
ment in sensitivity. On the other hand, if the noise corre- 
lation time is much shorter than the dephasing time T£ , 
then l/TJ is proportional to the square of the noise am- 
plitude [El], and thus T 2 *({0, 1}) = 4T 2 *({1, -1}). Con- 
sequently, there is no sensitivity improvement using the 
{1,-1} basis in a rapidly fluctuating environment. 



Photon-shot-noise limited sensitivity 

Due to finite photon collection efficiency and imper- 
fect spin-state readout, we do not actually realize the 
spin-projection-noise limited sensitivity given in Eq. Q. 
Instead, the experimental sensitivity is limited by the 
photon shot noise. We estimate the sensitivity by keep- 
ing the free precession interval, r w T 2 /2, constant and 
varying the applied longitudinal magnetic field, B z . The 
calibration of B z as a function of the current applied to 
our Hclmholtz coil is determined by varying the current 
and monitoring peaks in the Ramsey spectra similar to 
the inset of Fig. 1(d). The resulting fluorescence signal 
is plotted in Fig. [Sllfor both {1,-1} and {0, 1} bases. 
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FIG. SI: Magnetometer calibration curves for the {0,1} 
basis (up) and {1,-1} basis (bottom) , respectively. 



The curves are well described by the function 



P(B Z 



T , 



a— [1 + i?sin (2i:Am s g e fiBB z T 

T 



(S9) 



where 9 is a fitted phase, 2R is the measurement con- 
trast, a is the mean number of photons collected during 
a single readout pulse, and T is the total measurement 
duration. The photon shot noise for a given measurement 
is y/P(B z ), and therefore the photon-shot-noise-limited 

sensitivity is r\ v h = j^pj^^VT. The sensitivity is max- 
imized for 2TrAm s g e ^BB z T + 9 = jtitt, where m is an 
integer, giving 



1 



Vph 



2-KAm a g e \iB^[r R\fa' 



(S10) 



Compared to Eq. Q of the main text, we find the rela- 
tion r] p h = Vmin / (R\/ a ) ■ Typical values for our magne- 
tometer operating near the ideal case of r w T 2 /2 (Fig. 
Sll are 2R w 0.25 and a w 0.01. 



Allan deviation 

To study the Allan deviation of the measured noise, 
we used the following formula, 



<t b {t) = V0.5(S T ((n + l)r) - B r (nr)) 2 , (Sll) 



3 



where r is the gate time, 7TT means average over all in- 
teger n, and 3 T (nr) is the mean value of the measured 
magnetic field for the gate time in the range nr < t < 
(n + 1)t. 
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